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Extension of the optical spectroscopic method of isotope analysis of gaseous mixtures 
of hydrogen and deuterium to mixtures of water and heavy water provides a direct method 
for determining the total water content of materials. A determination can be made on a 
sample of less than 0.1 milliliter in less than 10 minutes to an accuracy of a few percent. 
Using a sample of at least 1 milliliter, an accuracy of the order of 0.1 percent is obtained in 
an hour or less. This paper describes the system and the investigation of optimum operating 
conditions and effect of variations in operating conditions on precision and accuracy of 
measurement, as well as limitations of the technique. A sample of blood or serum can be 
measured directly almost as easily as one of distilled water. The optical spectroscopic 
method can be used to determine the water content of most materials, whether in the solid, 
liquid, or gas phase. 



1. Introduction 

Following initial work on an optical spectroscopic 
analysis of gaseous mixtures of hydrogen and deu- 
terium [1,2], 2 a similar technique was tried with 
mixtures of water and heavy water. A preliminary 
study [3] indicated that it was possible to carry out 
such analyses to an accuracy of better than 3 percent. 
The work discussed in the present paper represents 
an extension of the preliminary study in an effort to 
improve both the precision and the range and to 
determine present limitations of the analysis. A 
detailed study was made to determine the influence 
of operating variables on precision and accuracy of 
measurement in order to find the optimum operating 
conditions. In these studies the following accuracy 
was obtained: (a) 0.1 percent of the intensity ratio, 
in the isotope ratio range 0.01 to 1.0; (b) 1 percent, 
in the range 0.00015 to 0.01 (see table 5). 

In addition to the above the rather practical prob- 
lems of developing the method into a routine analyti- 
cal technique were taken into consideration. For 
this reason efforts were made to achieve a rapid anal- 
ysis, and operating conditions were viewed with the 
possibility of constructing a simple flow system and 
optical system for routine use. Several types of 
routine measurements were made concurrently with 
the study of the technique. 

The technique used is as follows: The vapor from 
an aqueous sample containing a mixture of H 2 and 
D 2 is pumped through a tube where a high-fre- 
quency elect rodelcss discharge is maintained. The 
light emitted from the discharge tube contains lines 
of the Balmer series of hydrogen and deuterium with 
a very small amount of background radiation. The 
emitted light is dispersed by a monochromator, and 
the relative intensities of the lines of the Balmer 
series of hydrogen and deuterium are recorded and 
measured. The ratio of the intensities of the two 
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lines is a monotonic function of the relative concentra- 
tion of the isotopes, and to a high degree of accuracy 
is equal to the ratio of the isotope concentrations. 
For precise work the measurement is a relative, 
rather than an absolute one, in that the relationship 
between intensity ratio and concentration ratio is 
obtained from calibration curves taken on stain hue I 
samples whose compositions are determined by inde- 
pendent measurements. 

The statistical terms [4, 5] used in this report are 
defined as follows: "Precision" refers to the variation 
encountered during readings taken in close succession 
under identical conditions. "Reproducibility" re- 
fers to the variation among readings when they are 
spread over a considerable time period. It is not a 
fixed quantity but depends somewhat on the time 
interval and the history of the apparatus in that 
interval. " Accuracy " refers to the variation of the 
measured concentration from the true value and is 
determined from the precision and reproducibility 
of the measurements and from the accuracy of the 
calibrated samples. The relationship between repro- 
ducibility and time determines how frequently it is 
necessary to take calibration points. These fac- 
tors are all dependent on the stability of operating 
conditions. 

The above variations are expressed either as stand- 
ard deviation or as the coefficient of variation. These 
terms apply either for single measurements or for 
averages. As the standard deviation of the average 
is dependent on the number of readings taken, this 
value is most useful when compared to other sets of 
the same number of readings. For the purposes of 
this study, all standard deviations of the average 
refer to sets of 10 readings. 

Stability of this system may be viewed in the fol- 
lowing manner. Consider the measured ratio of 
intensities of the DjSBalmer line to the H/ffialmer 
line at 4860.0 A and 4861.3 A, respectively, and rep- 
resent this ratio by Di/Hi. This ratio is a function 
of a large number of variables such as the actual 
concentration ratio in the sample D/H, the tempera- 
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ture, the pressure, the flow rate, etc., and may be 
represented in the following fashion: 



D,/H 1 =D/H /(D/H, X u X, 



X n ), 



where Xi . . . are independent variables. In 
practice, operating conditions are found such that 
/ is very nearly unity. 

A completely stable set of operating conditions is 
one in which Di/Hi is a function of D/H only, and 
is unaffected by small changes of X x , X 2 . . . X n . 
From a practical point of view, the conditions of 
stability are that changes of Di/Hi due to variations 
of the X'$ are small compared with changes due to 
variations of D/H. Many of the experimental 
curves shown in the following sections are graphs of 
the ratio Di/Hi as a function of one of the independ- 
ent variables, with all other variables kept as constant 
as possible. 

The samples employed in these experiments are 
in all cases aqueous solutions or suspensions in which 
the solvent is a mixture of hydrogen oxide, H 2 0, and 
deuterium oxide, D 2 0. The number used to describe 
concentrations is, unless otherwise specified, the 
atom ratio of deuterium to hydrogen, D/H. This 
figure depends only upon the number of atoms of 
hydrogen and deuterium in the sample and is inde- 
pendent of the equilibrium between H 2 0, HDO, and 
D 2 0. 

In taking measurements of isotope ratios the 
vapor is pumped from the liquid through the dis- 
charge tube. This procedure raises two questions; 
what is the relationship between D/H in the liquid 
and the vapor and what fractionation (isotope sepa- 
ration) occurs in the flow process? The first ques- 
tion is purely academic with respect to the analysis 
of water. A given ratio of D to H in the liquid will 
always give rise to a reproducible ratio in the vapor, 
and as all measurements are taken relative to cali- 
bration samples, the measurement taken is, in all 
cases, referred back to the liquid. 

Because of the difference in vapor pressure between 
H 2 0, HDO, and D 2 0, the process of evaporating 
liquid from samples involves some degree of frac- 
tional distillation. This has been investigated experi- 
mentally to study the over-all effect and will be 
discussed in a later section. 

The processes occurring in the breakdown of water 
vapor in the radio-frequency field are imperfectly 
understood and will be discussed only briefly. The 
discharge through the water vapor in the tube 
causes excitation, ionization, and dissociation, giving 
rise to the following atoms and molecules from a 
mixture of normal and heavy water: 

H, D, HD, H 2 , D 2j OH, OD, H 2 2 , D 2 2 , H 2 0, HDO, 
D 2 0, HD0 2 , 

as well as other groupings, most of which may exist 
as normal, excited, or ionized molecules. The emis- 
sion spectrum from the discharge contains lines and 
bands from many of these components. The Balmer 
lines of hydrogen and deuterium are prominent in 
this spectrum. 



The p line of the Balmer series was chosen because 
the largest output signal was obtained for this line. 
This is a result of the product of the actual intensity 
and the phototube sensitivity being a maximum in 
the region of the ft Balmer line of the hydrogen 
spectrum. The a lines at 6561.0 A and 6562.8 A 
for D and H, respectively, also can be used; this pair 
has the advantage of a somewhat larger separation 
than the ft pair. (1.8 A compared to 1.3 A). 

2. Equipment 

2.1. Glass Flow System 

A glass flow system was employed to provide a 
steady flow of sample vapor through the discharge 
tube. A flow system decreases problems of contami- 
nation by previous samples and simplifies the prepa- 
ration of the sample for analysis. The system, 
shown photographically in figure 1 and schematically 
in figure 2, consists of the following components: 

1 . Sample bottles. The sample bottles are equipped 
with an inverted standard taper ground-glass joint 
so that rapid change of samples can be made. 

2. Sample connectors. The connectors are stand- 
ard taper ground-glass joints corresponding to those 
on the bottles. A special vacuum line (fig. 2) pro- 
vides a means of evacuating the air from the sample 
chamber before admitting sample into the main flow 
system. 

3. Pressure regulator. In order to vary and con- 
trol pressure in the discharge tube a parallel network 
is provided, consisting of a capillary and a continu- 
ously variable metal leak valve. In most of the 
studies described in this paper, the variable metal 
leak valve was used because fine pressure control 
and frequent adjustments were necessary. How- 
ever, for routine measurements a fixed nonmetallic 
orifice is preferable. 

4. Pressure gage. To measure pressure in the 
system, a pressure line connects through stopcocks 
to various critical locations. Pressure measurement 
is made by means of a commercial ionization gage, 
which utilizes alpha particles from a radium source 
to provide ionization. The gage can be used from 
pressures of 1/z to 1 cm. Pressures used in this report 
are gage readings. Absolute pressure of water vapor 
is found by multiplying by 1.163, the factor supplied 
by the gage manufacturer. 

5. Discharge tube. The discharge tube consists of 
a glass inner cylinder surrounded by a glass water 
jacket. Water flows through the jacket at a constant 
rate controlled by a regulator at the tap. Air 
bubbles in the water, which might affect intensity 
measurements, are eliminated by an extractor shown 
at point G in figure 1. 

6. Pumps. The two main vacuum lines shown 
in figure 2 connect with independent pumping 
systems. The left-hand line goes to a trap, single 
stage oil diffusion pump and forepump while the 
right-hand line goes only to a trap and mechanical 
pump. The diffusion pump was used in these 
investigations because a wide variation in conditions 
was to be studied. 
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Figure 1. Flow system for H— D analysis. 

A, Sample bottles; B, sample connectors; 0, metal leak valve; D, glass capillary; 
E, pressure sage; F, discharge tube; G, bubble extractor; H. trap; J, shield for 
high-frequency «>urce; K, lens; L, chopper; M, photomultipliejr tube. 



2.2. High-Frequency Source 

The energy necessary to excite the spectrum of 
the sample was provided by a commercial 150-Mc 
150-w oscillator. The energy is carried through a 
radio-frequency cable, which terminates in a tunable 
quarter-wavelength tank circuit with antennas held 
near the discharge tube. This terminal circuit is 
shown diagramatically in figure 2. A radio-frequency 
shield (point J in fig. 1) cuts down stray radiation, 
which might affect other electronic equipment. 

5Tyj • 2.3. Monochromator 

r 

The apparatus used to analyze the spectrum con- 
sists of a high-resolution grating monochromator 
with a photomultiplier detector and associated 
amplifier and recorder. 3 A detailed discussion of the 
original monochromator is given by Fas tie [6]. 

The plane diffraction grating used in the mono- 
chromator is 3 by 3 inches and is ruled with 30,000 
lines per inch; the instrument has a practical resolv- 
ing power (X/AX) of 55,000. This resolution is 
sufficient to observe some of the fine structure split- 
ting of H a and D a lines. The choice of optimum 
resolving power is governed by two considerations: 
it must be sufficient to resolve the lines to be analyzed 

3 This equipment was on field trial loan to the National Bureau of Standards 
_jom the Research Division of the Leeds & Nor thru p Co. 
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Figure 2. Schematic diagram of flow system. 
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Figure 3. Method of measuring intensities (Hj, D„ Gj). 

Height of H line (H is approximately 500 times that of H ghost line (G j) . The 
zero line is arbitrarily taken 6.6 A from the H line, and intensities are measured as 
distances above the zero line. From left to right, traces are of 0.4-percent and 10- 
percent D2O samples, respectively. 



and at the same time be low enough to allow for the 
integrated intensity of a line to be approximated by 
the peak height of that line. The resolving power 
was varied by varying the width of the entrance and 
exit slits until the p line of deuterium could be 
resolved in a normal water discharge (D/H= 0.000 15) 
[7]. Figure 3 shows the H, D, and H ghost lines 
resolved, and indicates, on the basis of the half 
width of the line, a resolving power of 20,000. This 
resolving power proved suitable for subsequent 
studies. 

The first Kowland ghost of the p line of hydrogen 
was spaced about 2.0 A from the line and was ap- 
proximately one five-hundredths as intense. The 
P line of deuterium appears between the H line and 
its first ghost. The ratio of ghost intensity to line 
intensity is a function of the optics of the system 
only, and under conditions of fixed optics may be 
taken as a constant. For low deuterium concentra- 
tions the ratio between the deuterium line and the H 
ghost is a fixed multiple of the ratio between the D 
and H lines. This fact proves useful in allowing 
comparison to be made between signals taken at 
the same scale of amplification, the calibration in 
this range being independent of the ratio of ampli- 
fication for different scales. 



2.4. Power Supply, Amplifier, and Recorder 

A chopper before the entrance slit of the mono- 
chromator converts the input signal to 60 cps. 
This signal is fed into a narrow band-pass syn- 
chronous amplifier, the output of which is recorded 
by a pen on 10-inch chart paper. The recorder 
has a response time of 1 second. The power supply 



for the photomultiplier is continuously variable from 
500 to 1,300 v. By using the photomultiplier with a 
sensitive linear a-c amplifier, which has an accurate 
amplification factor ranging from 2 to 800, very 
small signals can be observed. 

3. Experimental Results 

The primary purpose of this study was to deter- 
mine the effects of numerous interrelated variables 
on the speed, accuracy, and sensitivity of the method. 
The variables include pressure in the flow system, 
flow rate, sample size, solute composition of the sam- 
ple, temperature, and oscillator characteristics. 
These variables will be considered in detail. The 
percentage change in ratio per unit change in variable 
was computed for each operating curve and can be 
used to determine how closely each variable must be 
controlled to achieve a desired precision (see table 7). 

The methods of measuring intensity ratios for the 
0- to 1-percent D 2 range and the 1- to 99-percent 
D 2 range are shown in figure 3. As previously dis- 
cussed, the deuterium to ghost ratio is used instead 
of the deuterium to hydrogen ratio in the 0- to 1- 
percent range. In the 1- to 99-percent range the D 
and H lines were measured directly at a lower amplifi- 
cation. In both concentration ranges the baseline, or 
background, is somewhat arbitrary and is taken at a 
fixed distance from the deuterium peak. This causes 
no error in the final measurement because the calibra- 
tion is made in the same manner as the measurement. 

3.1. Contamination by a Previous Sample 

The speed with which a measurement can be taken 
is largely dependent on the time required to clean 
the system of one sample and saturate it with another 
so that a constant ratio is obtained. Metal parts 
are troublesome, but the metal leak valve usually 
was employed to gain fine pressure control in spite 
of its tendency to retain water with which it has 
previously been in contact. 

Figure 4 shows a series of experiments under ex- 
treme conditions in which a pure water sample was 
run and then switched to commercial heavy water. 
Curve a shows the effect, under the slowest condi- 
tions, of using the metal leak valve with no interme- 
diate pumping or flushing. Curve b shows the effect 
of using a glass capillary tube and a glass stopcock 
for pressure control with no intermediate pumping 
or flushing. For routine measurement the capillary 
tube and stopcock were used under the following con- 
ditions. The system was pumped with and without 
sample between extreme pressures in the ratio of at 
least 15 to 1. These changes were made rapidly for 
approximately 2 min, and then the pressure was 
stabilized to 400^- The effect of contamination by 
the previous sample was thereby reduced to that 
shown in curve c. 

These curves show that an all-glass system, with 
intermediate flushing, most quickly reduces contami- 
nation. Even with an extreme change in concentra- 
tion, equilibrium is reached to within 1 percent in 
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Figure 4. Residual contamination. 

Pressure equals 400 ». a, Metal leak valve, with no intermediate pump- 
ing or flushing between samples; b, glass stopcock, with no intermediate 
pumping or flushing; c, glass stopcock, with intermediate pumping and 
flushing. 

less than 3 min and complete equilibrium in approxi- 
mately 1 ln\ The flow system shown in figure 2 was 
designed for maximum variability, but for routine 
use, a much simpler system is advantageous. Also 
the elimination of side projections should reduce the 
time to remove the previous sample. 

3.2. Pressure 

The operating pressure in the discharge tube has 
an effect on the intensity ratio, especially in the 0- to 
1-percent D 2 range (see fig. 5). Flow rate was not 
held constant, so actually the traces show ratio versus 
a combination of pressure and flow rate at constant 
pumping speed as determined by settings of the 
metal leak valve and the speed of the pump. The 
flow rate in milligrams per minute decreased from 
0.011 at a pressure of 900m to 0.0023 at 400/x and 
0.001 at 200/i. The pressure gage used had the dis- 
advantage of large metal surfaces, resulting in appre- 
ciable delay in removing traces of previous samples, 
an effect that is particularly undesirable when the 
concentration is being changed frequently. There- 
fore, the gage was isolated from the system after 
initial pressure setting and during the course of most 
intensity-ratio measurements. The disadvantage of 
not taking continuous pressure readings was over- 
come by using the discharge length as a pressure 
indicator because it is a sensitive function thereof. 
Although this variation could not be used for absolute 
measurement of pressure, it provided a good check 
on the constancy and reproducibility of pressure and 
flow conditions. 

Intensity as a function of pressure was investigated 
at various concentrations, and the results are shown 
in figure 6. The curves are quite similar at all con- 
centrations. An operating pressure of approximately 
400 /x was chosen for routine use because it makes use 
of maximum intensity and the flattest portion of the 
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Figure 5. Intensity ratio versus pressure — seven 
samples at various concentrations. 

The calibrated concentration of each sample by weighing is 
given at end of each trace. 
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Figure 6. Total H^-r-D^ intensity versus pressure at various 
concentrations. 



intensity versus pressure curve. It also occurs at 
regions of the ratio versus pressure curves that are 
reasonably flat. The percentage change of ratio per 
unit change in pressure at 400^ as computed from 
figure 5 is 0.020 percent per micron. This means 
that the pressure must be held to ±5u to achieve a 
precision of ±0.1 percent. 
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3.3. Fractionation 

At 20° C there is a difference of 2.3 mm of mercury 
between the vapor pressure of H 2 and D 2 0, with 
HDO having an intermediate pressure [7]. Thus a 
sample that is being evaporated will begin to concen- 
trate deuterium. To check on the error to be ex- 
pected from this effect, a 0.5-ml sample was run 
continuously to dryness, with frequent checks made 
on the D/H ratio. The results at an operating 
pressure of 200 p. are shown in figure 7. From this 
curve the minimum sample size can be computed for 
a given flow rate and allowable error. This is shown 
in detail in section 5. 
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Figure 7. Fractionation effect. 
Flow was continuous at 0.063 ml/hr. Di/Gi scale is greatly expanded. 

3.4. Temperature 

The ratio of the vapor pressure of H 2 0, D 2 0, and 
HDO changes with temperature, so that some error 
might be expected from variation of sample tempera- 
ture. To study this effect, sample bottles were 
immersed in a constant-temperature bath during 
these determinations. Table 1 shows the results 
obtained at temperatures ranging from 0.4° to 25.8° 
C. From this table the percentage change in ratio 
per unit change in sample temperature was com- 
puted to be 0.15 percent per deg C. 

Table 1. Intensity ratio versus sample temperature 
Pressure equals 500 (i. 



Sample 


Tempera- 
ture 


Ratio 
Di/Hi 


S-4. 


C 

J 4.0 
I 25.8 
( 0.4 
3.0 
1 17.0 
I 22.0 


0. 3895 
.3965 
.0454 
. 0452 
.0466 
.0470 


S-9 





3.5. Flow-Rate 

Flow rate affects the D/H ratio, the rate of frac- 
tionation, and the size of the smallest sample 
measurable. It is a function of pressure, the 
geometry of the flow system, and the speed of the 
pump. Table 2 gives data obtained for Di/Hi ratio 



as a function of the flow rate at constant pressure. 
This table shows that Di/Hi decreases with de- 
creased flow rate and indicates that at least part of 
the effect of change of pressure (fig. 5) is caused by 
the reduced flow rate at lower pressures. The 
percentage change of ratio per unit change of flow 
rate computed from table 2 is 0.5 percent per gram 
a minute. 

Table 2. Intensity ratio versus flow rate 

Pressure equals 400 m- Flow rate obtained by weighing sample before and after 

flow. 



Flow 


Ratio 
Di/Hi 


g/min 

0.01125 

. 00983 

. 00331 


0.899 
.894 

.857 



3.6. Oscillator and Discharge Characteristics 

The discharge characteristics can be varied by 
changes in oscillator plate current, antenna tuning, 
position of electrodes with respect to the tube, and 
the temperature of the cooling stream in the jacket 
of the discharge tube. Plate current is variable 
between 100 and 200 ma, and the tuning is adjusted 
by the tuning capacitor in figure 2 as well as by the 
position of the radio-frequency shield and electrodes. 
The radio-frequency shield and electrodes were placed 
to obtain maximum intensity, and the tuning capaci- 
tor was varied to give a range of discharge lengths. 
Ratio of Di/Hi as a function of discharge length is 
shown in table 3. The effect on Di/Hi of changing 
the oscillator plate current is shown in table 4. It 
is seen that there is no appreciable effect of the 
discharge length or plate current on the measured 
ratio. The temperature of the cooling stream was 
varied by changing the water flow rate through the 
discharge-tube jacket. Di/Hi varied less than 1 
percent over the extreme range of temperature from 
a stream of the coldest water available from the 
tap (approximately 18° C) to no cooling (greater 
than 200° C). 

Table 3. Intensity ratio versus percentage maximum discharge 
length 

Pressure equals 500 n 



Sample 


Percentage of 
maximum dis- 
charge length 


Ratio 
Di/Hi 


S-9 


f 100 
J 88.6 
1 73.8 
I 49.4 
( 100 
98 

91.8 
/ 82.7 
75.1 
64.7 
52.0 
49.0 


4.71 

4.78 
4.85 
4.48 
49.74 
49.5 
49.71 
49.79 
49.97 
50.17 
50.04 
49.79 


S-12 





Table 4. Intensity ratio versus oscillator plate current 
Pressure equals 500 n 



Sample 


Plate 

current 


Ratio 
Di/Hi 


S-9 


via 

100 

122 

140 

175 

. 190 

( 100 

125 

130 

150 

185 

200 

I 210 


0. 0503 
.0499 
.0491 
.0491 
.0490 
1.015 
1.016 
1.013 
1.015 
1.012 
1.012 
1.012 


S-12 





3.7. Solute Composition 

As the optical spectroscopic method of analysis 
may be extended for use with many aqueous solu- 
tions, preliminary investigation was made of the 
variation of Dt/Hi ratio with the amount of solids 
dissolved in the sample. This test was made with 
NaCl. With concentrations of percent, 1 percent, 
10 percent, and a saturated solution, no variation of 
Di/Hi was found within the experimental error. 



4. Calibration Curves, Precision, and Repro- 
ducibility 

Calibration curves were necessary to achieve the 
highest accuracy. Figure 8 shows the 0- to 1 -percent 
D 2 range and figure 9 the 1- to 80-percent D 2 
range. 4 Samples were made bv mixing portions of 
distilled water and commercial D 2 (99.8%) that had 
been weighed to an accuracy of 0.1 mg. 

Table 5 gives examples of the isotope content ex- 
pressed, for comparison, in ratio and in percent, 
and shows the precision obtained under optimum 
conditions. 



Table 5. Comparison of ratio method and percentage concen- 
tration method of expressing isotope content 





Standard 


Coefficient 


Deuterium 


Standard 
deviation of 
average of 
concentra- 
tion » 


Coefficient 


Ratio 


deviation 


of 


concentra- 


of 


D/H 


of average 


variation 


tion 


variation of 




of ratio a 


of average b 


100D/D+H 


average b 








% 






500 


5.0 


1.0 


99.8 


0.002 


0.002 


100 


0.1 


0.1 


99.0 


.001 


.001 


10 


ul 


.1 


90.9 


.01 


.01 


1 


.001 


.1 


50.0 


.25 


.05 


.1 


. 0001 


.1 


9.09 


.01 


.01 


.01 


.00001 


.1 


0.99 


.001 


.1 


.001 


.000007 


.7 


.100 


.0007 


.7 


.0001 


. 000002 


2.0 


.01 


.0002 


2.0 



a Based on a Se1 of 10 readings. 

b Standard deviation of average, expressed in percent. 



4 Calibra 1 ton curves were taken on samples up to 99.8 percent of D2O (commer- 
cially available), but do extensive measurements were made in the high-concen- 
tration range. The analysis in this range differs from that in the low ranges 
only in the method of preparing standard samples. 
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Figure 8. Calibration curve for 0- to 1-percent D 2 
concentration range. 

Samples were made up by weighing. 



D 2 CONCENTRATION , percent 
33 50 67 75 



3.0 



2.0 



x 




D 2 0/H 2 (WEIGHED) 

Figure 9. Calibration curve for 1- to 80-percent D 2 
concentration range. 

80- to 99-percent portion of the curve is omitted because of the extended scale. 
Samples were made up by weighing. 



Reproducibility was measured on a 21 -percent 
D 2 sample and a 0.2-percent D 2 sample. Meas- 
urements were taken at varying time intervals, and 
in some cases with major changes taking place 
during the interval, both in the apparatus and in 
the concentration ratios. Results are shown in 
figure 10 for 4 portions of the 21-percent sample. 
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S-29 
SAMPLE NUMBER 



Figure 10. Precision and reproducibility for 21-percent 
D 2 samples. 



The data were taken over a period of 3 months, and 
dates are given below, reading from left to right 
on figure 10. 



Sample 


Date 


Sample 


Date 




[Dec. 


9, 1952 




[Dec. 4, 1952 


S-28 


Dec. 


9, 1952 


S-30 


Dec. 4,1952 




Dec. 


23, 1952 




[Dec. 23, 1952 




Dec. 


1, 1952 


S-31 


/Dec. 19, 1952 
IDec. 19, 1952 




Dec. 


1, 1952 




Dec. 


2, 1952 








Dec. 


2, 1952 








Dec. 


3, 1952 








Dec. 


3, 1952 






S-29 


Dec. 
Dec. 
Dec. 
Dec. 
Dec. 
Dec. 
Feb. 


4, 1952 

5, 1952 
5, 1952 
5, 1952 
5, 1952 

22, 1952 
25, 1953 







Each point shows the average of 10 readings, 
with the precision depicted as vertical spread. The 
grand average is a solid horizontal line, and the 
standard deviation computed directly from all 21 
averages is shown by dotted lines on either side of 
the grand average. Thus the relation between 
precision and reproducibility at this stage of develop- 
ment can be seen at a glance. The results obtained 
for the 0.2-percent sample were similar. Data for 
it are listed in table 6. 



Table 6. Reproducibility of intensity ratio for a 0.2-percent 
D 2 sample over a 1-month period. 

Each ratio is the average of 10 readings. Pressure equals 200/z. 



Date, 1953 


Di/Gi a 


Date, 1953 


Di/Gi a 


April 28 


0.760 
.737 
.736 
.733 
.729 
.724 
.744 
.751 
.743 


May 11 


.753 
.740 
.769 
.754 
.765 
.753 


May 1 


May 15 

May 18 


May 2 


May 4 . ... 


Do 

May 19 


May 7 

Do 

Do 


May 20. _ 


Average... ._. . _. 


May 8 


0.746 


Do 







5. Discussion 

The great number of variables that influence the 
measurement of D/H ratios makes it impossible to 
be certain that the final operating conditions are in 
fact the optimum ones. Indeed, it is not necessary 
to assume that there is a unique set of optimum 
conditions for this analysis. The operating condi- 
tions chosen therefore are the "best set", based upon 
data available and the nature of the measurements 
to be performed. Some of the final data used to 
evaluate operating conditions are shown in table 7. 
This table illustrates that there is no extreme 
problem in keeping the operating conditions constant 
enough so that the measurement can be made within 
the desired accuracy. 

Table 7. Operating variables 



Operating condition 
varied 


Percentage change of 
ratio per unit change 
of variable 


Routine operating 
conditions 


Sample temperature- 
Pressure 

Salt concentration--. 
Flow rate. . 


0.15 percent per deg C._. 
0.020 percent per micron . 
percent per mole/liter .. 
0.5 percent per g/min.. 

percent per percentage 
change in discharge 
length. 

percent per ma 

0.005 percent per deg C„. 


Room temperature. 

400 microns. 

Any. 

Depending on sample 
size and desired 
speed of measure- 
ment. 

Maximum discharge 
length. 

175 ma. 

Room temperature. 


Tuning- 


Oscillator plate cur- 
rent. 

Discharge tube tem- 
perature. 



» Gj is the measured intensity of the ghost of the H line. 



Accessory factors enter into the choice of best 
conditions. The fractionation curve (fig. 7) can. be 
represented as RS k =K, where R is the isotope ratio, 
S is the sample size, and k and K are experimental 
constants. In differential form the fractionation 
curve can be represented as dR/R= — (k/S)dS, where 
the left-hand term represents the percentage change 
of ratio due to fractionation. Setting dR/R equal 
to 6, the maximum tolerable fractionation error, 
and realizing that dS=—fAt, the flow-rate times 
the time increment, we arrive at the following ex- 
pression for minimum sample size: 



S= 



omsfAt 



The value of k, 0.068, was obtained from the data 
shown in figure 7. As has been mentioned previ- 
ously, flow-rate depends on operating pressure, flow- 
geometry, and pumping speed, and the time of 
measurement depends on the time required to clean 
the system of the previous sample and to equalize 
the pressure. Hence, sample size is some complex 
function of the major operating variables. Under 
normal conditions at a 400-^t operating pressure 
the flow rate was 0.0023 ml/min, and the time per 
determination was 20 min. The minimum sample 
size is therefore 0.003 IS/e, or 0.3 ml, for a 1 -percent 
error. The absolute minimum sample size under 
present conditions involves an operating pressure 
of 200 }i and a maximum error of 3 percent. • Under 
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these conditions, a 7-min measurement is possible, 
and the minim um sample size is 0.02 ml. 

The operating conditions that were found to be 
most suitable with the equipment described are 
shown in table 7. 

This study is not exhaustive, inasmuch as there 
are many combinations of variables that were not 
tried. In addition, no attempt was made to vary 
systematically such things as the oscillator frequency 
or discharge-tube geometry. A few preliminary 
experiments with a 2,240-Mc oscillator indicate that 
a stable set of operating conditions can be set up at 
that frequency. 

As pointed out in the beginning of this paper, the 
chief object of these studies was to determine how 
to minimize error, sample size, and time of determi- 
nation. Of course, it is impossible to minimize all 
three factors at once. Highest precision requires 
many determinations per sample, hence a long time 
and large sample size. The minimum values attain- 
able at the present stage of development and the 
maximum values obtained while minimizing the 
others are given as follows: 



Factor 


Minimum 


Maximum 


Precision (coefficient of 
variation of the aver- 
age) _ _ 

Sample size 


0. 1 to 2 per- 
cent. 
0. 02 ml 
7 min 


3 percent 

1 ml 
60 min 


Time per sample 



With respect to precision and time per determina- 
tion, the optical spectroscopic case compares favor- 
ably with other methods. It finds particular appli- 
cation where the water samples contain impurities. 
Whereas other methods require one or two stages 
of purification, the sample vapor can be analyzed 
directly by this method. This method appears to be 
considerably more rapid and, for similar accuracy, 
less expensive than oilier methods. 

While it is not the function of this paper to discuss, 
in detail, applications of hydrogen-deuterium analy- 
sis, we shall review briefly in the following para- 
graphs some of the most important uses to which the 
optical spectroscopic method has been placed. The 
dilution technique for (let ermining the total water 
content of materials involved in these uses is de- 
scribed in reference [3]. Essentially, it consists in 
dissolving or suspending a known amount of the 
material to be analyzed in a mixture of D 2 and 
H 2 and then determining spectroscopically the 
resulting change in the ratio of D 2 to H 2 0. 

1 . Determination of D 2 content of blood. In this 
application the blood, serum, or plasma is placed in 
the sample bottle and the vapor analyzed in the 
usual manner. No prior preparation of the blood is 
necessary and high precision is obtainable. This 
technique finds application in the determination of 
total body water by D 2 dilution and in the study of 
water exchange in capillary beds. Over 150 blood 
samples have been run, and no difficulty has been 



encountered in adapting the techniques to routine 
analysis. 

2. Determination of water of hydration or water of 
crystallization. The crystal to be analyzed is 
dissolved in D 2 0, and the final intensity ratio gives 
a measure of the H 2 originally in the crystal. 
Several hydrates have been tried, and the results 
agree with other methods. 

3. Determination of labile hydrogen. A weighed 
quantity 7 of dry material containing labile hydro- 
gen is dissolved in D 2 and the resultant ratio 
measured . From the number of exchangeable hydro- 
gen atoms per molecule, information about structure 
of the molecule is obtained. 

The most extensive work has been done with blood, 
for which a brief description of the 4 routine applica- 
tion follows. 

Army Medical Service gave injections of 100 g of 
commercially pure D 2 to several hemorrhagic fever 
patients and sampled blood at various time intervals. 
The serum samples were analyzed at the National 
Bureau of Standards to determine the total body 
water content of the patients. Seventy-three sam- 
ples in the concentration range 0.029 to 0.425 percent 
gave precisions of the order of 0.7-percent coefficient 
of variation of the average. Ten readings were 
taken per sample, using the metal leak valve for 
pressure control. The time per sample ranged 
between 25 and 55 min, depending on stability of 
operating conditions. 

Reference [8] gives details of a study of trans- 
capillary exchange rates of D 2 and thiocyanate in 
the forearm of man. A study has since been made of 
the exchange rates from the artery of one arm to the 
vein of the other. In this case only five readings 
were taken per sample and metal parts were 
eliminated, thus cutting the total time per sample 
down to 15 min, plus data calculation time. Th< 
precision of measurement was of the order of 1.5- 
percent coefficient of variation of the average for 
17 samples in the concentration range 0.015 to 2.00 
percent of D 2 0. 
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